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Abstract 

Rotating waves are periodic solutions in 5*0(2) equivariant dynamical systems. Their pre- 
cession frequency changes with parameters and it may change sign, passing through zero. When 
this happens, the dynamical system is very sensitive to imperfections that break the SO (2) 
symmetry and the waves may become trapped by the imperfections, resulting in steady solu- 
tions that exist in a finite region in parameter space. This is the so-called pinning phenomenon. 
In this study, we analyze the breaking of the 5*0(2) symmetry in a dynamical system close to a 
Hopf bifurcation whose frequency changes sign along a curve in parameter space. The problem 
is very complex, as it involves the complete unfolding of high codimension. A detailed analysis 
of different types of imperfections indicates that a pinning region surrounded by infinite-period 
bifurcation curves appears in all cases. Complex bifurcational processes, strongly dependent on 
the specifics of the symmetry breaking, appear very close to the intersection of the Hopf bifur- 
cation and the pinning region. Scaling laws of the pinning region width, and partial breaking 
of 50(2) to Z m , are also considered. Previous and new experimental and numerical studies of 
pinned rotating waves are reviewed in light of the new theoretical results. 



1 Introduction 

Dynamical systems theory plays an important role in many areas of mathematics and physics 
because it provides the building blocks that allow us to understand the changes many physical 
systems experience in their dynamics when parameters are varied. These building blocks are the 
generic bifurcations (saddle- node, Hopf, etc.) that any arbitrary physical system experiences under 
parameter variation, regardless of the physical mechanisms underlying the dynamics. When one 
single parameter of the system under consideration is varied, codimcnsion-one bifurcations are 
expected. If the system depends on more parameters, higher codimension bifurcations appear and 
they act as organizing centers of the dynamics. 

The presence of symmetries changes the nature and type of bifurcations that a dynamical 
system may undergo. Symmetries play an important role in many idealized situations, where 
simplifying assumptions and the consideration of simple geometries result in dynamical systems 
equivariant under a certain symmetry group. Bifurcations with symmetry have been widely studied 
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symmetries are only approximately fulfilled, and the breaking of the symmetries, due to the presence 
of noise, imperfections and/or other phenomena, is always present. There are numerous studies of 
how imperfect symmetries lead to dynamics that are unexpected in the symmetric problem, e.g. 
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). However, a complete theory is currently unavailable. 



One observed consequence of imperfections in systems that support propagating waves is that 



the waves may become trapped by the imperfections (e.g., see Keener 1987 Westerburg fc Busse 
2003 Thiele & Knobloch 2006a|& ). In these various examples, the propagation direction is typically 



biased. However, a more recent problem has considered a case where a rotating wave whose sense of 
precession changes sign is pinned by symmetry-breaking imperfections (Abshagen et al. 2008). We 
are unaware of any systematic analysis of the associated normal form dynamics for such a problem 
and this motivates the present study. 

When a system is invariant to rotations about an axis (invariance under the SO(2) symmetry 
group), 50(2)-symmetry-breaking Hopf bifurcations result in rotating waves, consisting of a pattern 
that rotates about the symmetry axis at a given precession frequency without changing shape. 
This frequency is parameter dependent, and in many problems, when parameters are varied, the 
precession frequency changes sign along a curve in parameter space. What has been observed in 
different systems is that in the presence of imperfections, the curve of zero frequency becomes a 
band of finite width in parameter space. Within this band, the rotating wave becomes a steady 
solution. This is the so-called pinning phenomenon. It can be understood as the attachment of the 
rotating pattern to some stationary imperfection of the system, so that the pattern becomes steady, 
as long as its frequency is small enough so that the imperfection is able to stop the rotation. This 
pinning phenomenon bears some resemblance to the frequency locking phenomena, although in the 
frequency locking case we are dealing with a system with two non-zero frequencies and their ratio 
becomes constant in a region of parameter space (a resonance horn), whereas here we are dealing 
with a single frequency crossing zero. 

In the present paper, we analyze the breaking of 5*0(2) symmetry in a dynamical system close 
to a Hopf bifurcation whose frequency changes sign along a curve in parameter space. The analysis 
shows that breaking 5*0(2) symmetry is much more complex than expected, resulting in a bifur- 
cation of high codimension (about nine). Although it is not possible to analyze in detail such a 
complex and high-codimension bifurcation, we present here the analysis of five different ways to 
break 50(2) symmetry. This is done by introducing into the normal form all the possible terms, 
up to and including second order, that break the symmetry, and analyzing each of these five terms 
separately. Three of these particular cases have already been analyzed in completely different con- 



texts unrelated to the pinning phenomenon (Gambaudo 1985 Wagener 2001 Broer et al.. 2008 



Saleh & Wagener 2010). In the present study, we extract the common features that are associated 



with the pinning. In all cases, we find that a band of pinning solutions appears around the zero 
frequency curve from the symmetric case, and that the band is delimited by curves of infinite-period 
bifurcations. The details of what happens when the infinite-period bifurcation curves approach the 
Hopf bifurcation curve are different in the five cases, and involve complicated dynamics with several 
codimension-two bifurcations occurring in a small region of parameter space as well as several global 
bifurcations. 

Interest in the present analysis is two-fold. First of all, although the details of the bifurcational 
process close to the zero-frequency Hopf point are very complicated and differ from case to case, for 
all cases we observe the appearance of a pinning band delimited by infinite-period bifurcations of 
homoclinic type that, away from the small region of complicated dynamics, are SNIC bifurcations 
(saddle-node on an invariant circle bifurcation, e.g. see Strogatz 1994). Secondly, some of the 



scenarios analyzed are important per se, as they correspond to the generic analysis of a partial 
breaking of the 50(2) symmetry, so that after the introduction of perturbations, the system still 
retains a discrete symmetry (the Zi case is analyzed in detail). 

The paper is organized as follows. In section ^2] the properties of a Hopf bifurcation with SO (2) 
symmetry with the precession frequency crossing through zero are summarized, and the general 
unfolding of the 50(2) symmetry breaking process is discussed. The next sections explore the 
particulars of breaking the symmetry at order zero (*j3]), one (<Q and two (Sj5|. Sections f^4]and { 5.1 
are particularly interesting because they consider the symmetry-breaking processes 50(2) — » Z2 
and 50(2) — > Z3 which are readily realized experimentally. Section ^extracts the general features 
of the pinning problem from the analysis of the specific cases carried out in the earlier sections. 
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Section SjSJpresents comparisons with experiments and numerical computations in two real problems 
in fluid dynamics, illustrating the application of the general theory developed in the present study. 
Finally, in <j9j conclusions and perspectives are presented. 



2 Hopf bifurcation with 50(2) symmetry and zero frequency 

The normal form for a Hopf bifurcation is 

z = z(fi + ito — c|z| 2 ), (1) 

where z is the complex amplitude of the bifurcating periodic solution, /i is the bifurcation parameter, 
and to and c are functions of fi and generically at the bifurcation point (/i = 0) both are different 
from zero. It is the non-zero character of to that allows one to eliminate the quadratic terms in z in 
the normal form. This is because the normal form z = P(z, z) satisfies (e.g., see Haragus &: Iooss[ 



2011) 



P{e~ lult z, e lult z) = e- Iwt P(z, z), (2) 

where P is a low order polynomial that captures the dynamics in a neighborhood of the bifurcation 
point. If to — 0, this equation becomes an identity and P cannot be simplified. The case to — is 
a complicated bifurcation and it depends on the details of the double-zero eigenvalue of the linear 
part L of P; as z — x + iy is complex, the matrix of L using the real coordinates (x,y) is a real 
2x2 matrix. If L is not completely degenerate, that is 



1 




(3) 



(4) 



then we have the well-studied Takens-Bogdanov bifurcation, whereas the completely degenerate 
case, 

'0 0^ 

v o o y 

is a high-codimension bifurcation that has not been completely analyzed. 
If the system has 150(2) symmetry, it must also satisfy 

P(e imf) z, e- irn0 z) = e im9 P{z 7 z), (5) 

where Z m is the discrete symmetry retained by the bifurcated solution. When the group Z m is 
generated by rotations of angle 27r/m about an axis of ra-fold symmetry, as is usually the case with 
SO(2), then the group is also called C m . Equations |2| and |5|) are completely equivalent and have 
the same implications for the normal form structure. Advancing in time is the same as rotating the 
solution by a certain angle {tot = m0); the bifurcated solution is a rotating wave. Therefore, if to 
becomes zero by varying a second parameter, we still have the same normal form ([!]), due to ([5]), 
with to replaced by a small parameter v. 

z — z{fi + iv — c|z| 2 ). (6) 

The Hopf bifurcation with SO{2) symmetry and zero frequency is, in this sense, trivial. Introducing 
the modulus and phase of the complex amplitude z = re 1 ^ , the normal form becomes 



r = r(/i — ar 2 ) 
(j) = v — br 2 , 



(7) 



where c = a + ib, and let us assume for the moment that a and b are positive. The bifurcation 
frequency in Q is now the small parameter v. The bifurcated solution RW m exists only for [i > 0, 
and has amplitude r — \fixfa and frequency to = v — bfi/a. The limit cycle RW m becomes an 
invariant set of steady solutions along the straight line fi — au/b (labeled L in figure [TJ) where 
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Figure 1: Hopf bifurcation with SO (2) symmetry and zero frequency; part (a) shows the bifurcation 
diagram, where the thick lines are bifurcation curves, and part (b) shows the bifurcations along the 
path A shown in (a). The fixed point curve is labeled with the signs of its eigenvalues. In regions 
II_ and II + the limit cycles, born at the Hopf bifurcations H_ and H + , rotate in opposite senses. 
L is the line where the limit cycle becomes an invariant curve of fixed points. 



the frequency of RW m goes to zero; the angle between L and the Hopf bifurcation curve (the 
horizontal axis fx = 0) is ao- The bifurcation diagram and a schematic of the bifurcations along 
a one-dimensional path is also shown in figure [T] The bifurcation point [i = v = 0, labeled ZF 
(zero-frequency Hopf point) in figure [lja), is a codimension-two bifurcation. It coincides with the 
generic Hopf bifurcation, except that it includes a line L along which the bifurcated solution has 
zero frequency. 

Assuming c 7^ 0, we can simplify |7]) by scaling z so that |c| = 1; we will write 



c = a + \b = ie 



sin ao + 1 cos ao, 



a"*0 



(8) 



which helps simplify subsequent expressions. The case a and b both positive, which we will analyze 
in detail in the following sections, corresponds to one of the fluid dynamics problems that motivated 

and 



the present analysis (see Abshagen et al. 2008 Pacheco et al. 2011 



i i.l). For other signs of 



a and 6, analogous conclusions can be drawn. It is of particular interest to consider the subcritical 
case a < as it corresponds to the other fluid dynamics problem analyzed here 



(sec 1 



Marques et al. 

2007 Lopez fc Marques] |2009 and i.2). By reversing time and changing the sign of \x and v, we 
obtain exactly the same normal form (171) but with the opposite sign of a and b. By changing the sign 
of <j> and v, we obtain ^ with the opposite sign of b. Therefore, all possible cases corresponding 
to different signs of a and b can be reduced to the case where a and b are both positive. 



2.1 Unfolding the Hopf bifurcation with zero frequency 

If the 50(2) symmetry in the normal form (|6| is completely broken, and no symmetry remains, 
then the restrictions imposed on the normal form by ^ disappear completely and all the terms in 
z and z missing from Q will reappear multiplied by small parameters. This means that the normal 
form will be 



i = z(fi + iu — c\z\ 2 ) + ei + + e$z 2 + e^zz + e^z 2 



(9) 



where additional cubic terms have been neglected because we assume c 7^ and that cz\z\ 2 will be 
dominant. As the Q are complex, we have a problem with 12 parameters. Additional simplifications 
can be made in order to obtain the so-called hypernormal form; this method is extensively used 



by Kuznetsov (2004), for example. Unfortunately, many of the simplifications rely on having some 



low-order term in the normal form being non-zero with a coefficient of order one. For example, if 
to =/= 0, it is possible to make c real by using a time re-parametrization. In our problem, all terms up 
to and including second order are zero or have a small coefficient, and so only a few simplifications 
are possible. These simplifications are an infinitesimal translation of z (two parameters), and an 
arbitrary shift in the phase of z (one parameter). Using these transformations the twelve parameters 
can be reduced to nine. In particular, one of either £4 or £5 can be taken as zero and the other can be 
made real. By rescaling z, we can make c of modulus one, as in (J8|. A complete analysis of a normal 
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form depending on nine parameters, i.e. a bifurcation of codimension of about nine, is completely 
beyond the scope of the present paper. In the literature, only codimension-one bifurcations have 
been completely analyzed. Most of the codimension-two bifurcations for ODE and maps have 



also been analyzed, except for a few bifurcations for maps that remain outstanding (Kuznetsov 



2004 1 . A few codimension-three and very few codimension- four bifurcations have also been analyzed 



(Chow et al. 1994 Dumortier et al. 1997), but to our knowledge, there is no systematic analysis 



of bifurcations of codimension greater than two. 

In the following sections, we consider the five cases, t\ to £5, separately. A combination of 
analytical and numerical tools allows for a detailed analysis of these bifurcations. We extract the 
common features of the different cases when <C \J /i 2 + v 2 , which captures the relevant behavior 
associated with weakly breaking 50(2) symmetry. In particular, the £2 case exhibits very interesting 
and rich dynamics that may be present in some practical cases when the SO(2) symmetry group is 
not completely broken and a Z2 symmetry group, generated by the half-turn 9 — > 8 + 7r, remains. 

Some general comments can be made here about these five cases, which are of the form 



z(p- 



(10) 



for integers < q < p < 2, excluding the case p = q = 1 which is 50(2) equivariant and so e can 
be absorbed into /j, and v. By changing the origin of the phase of z, we can modify the phase of e 
so that it becomes real and positive. Then, by re-scaling z, time t, and the parameters /j, and v as 

in) 



(z, t, n, v) -> (e d z, e~ 25 t, e M fi, e 2 V), S = 



2(5. 



,25, 



P 



we obtain (10) with e = 1, effectively leading to codimension-two bifurcations in each of the five 
cases. We expect complex behavior for /i 2 + iP < e 2 , when the three parameters are of comparable 
size, while the effects of small imperfections breaking 50(2) will correspond to /i 2 + v 2 ^S> e 2 . 
From now on e = 1 will be assumed, and we can restore the explicit e-dependence by reversing 



the transformation (11). Three of the five normal forms (11) have been analyzed in the literature 



(discussed below), focusing on the regions where [A, v and e are of comparable size; here we will also 
consider what happens for /i 2 +^ 2 3> £ 2 which is particularly important for the pinning phenomenon. 

The normal forms corresponding to the e%, £2 and £3 cases have already been analyzed in con- 
texts completely different to the 50(2) symmetry-breaking context considered here. The context 
in which these problems were studied stems from low-order resonances in perturbed Hopf problems. 



Gambaudo (1985) studied time-periodic forcing near a Hopf bifurcation point, analyzing the prob- 



lem using the Poincare stroboscopic map. The normal forms corresponding to the 1:1, 1:2 and 1:3 
strong resonances coincide with the normal forms we present below for cases with only the £1, £2 
and £3 terms retained in ([9]), respectively. Later, motivated by a problem of a nonlinear oscillator 
with damping and quasi-periodic driving, a series of papers extended the strong resonances results 



of Gambaudo ( 1985 ) by studying the semi-global bifurcations for periodically and quasi-periodically 



perturbed driven damped oscillators near a Hopf bifurcation (see Wagener 2001 Broer et al. 2008 



Saleh & Wagener 2010 and references therein). The other two cases we consider, with only the £4 
or the £5 terms retained in ^ , do not appear to have been studied previously. They fall outside of 
the context in which the other three were studied because they do not correspond to any canonical 
resonance problem. We should point out that within the resonance context, the three cases studied 
would not make sense to consider in combination (they correspond to completely distinct frequency 
ratios and so would not generically occur in a single problem). In contrast, within the context 
motivating our study, all five cases correspond to different ways in which the 50(2) symmetry of 
a system may be broken, and in a physical realization, all five could co-exist. In the following 
sections, we present a detailed analysis of all five cases. 



3 Symmetry breaking of SO (2) with an e term 



The normal form to be analyzed is ( 10 ) with p = q = and e = 1: 

i = z(n + iv — c|z| 2 ) + 1. 



(12) 
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Figure 2: (a) Coordinates (u,v) in parameter space adapted to the line L, which coincides with 
zero frequency curve in the unperturbed SO(2) symmetric case, (b) Steady bifurcations of the fixed 
points corresponding to the normal form ( |12[ ). SN± and SNo are saddle- node bifurcation curves 
and Cusp-t are cusp bifurcation points. In region III there exist three fixed points, and only one in 
the rest of parameter space. 



This case has been analyzed in Gambaudo (1985); Wagener (2001); Broer et al. (20081; Saleh & 



Wagener (2010) 



It is convenient to introduce coordinates (u, v) in parameter space, rotated an angle ao with 
respect to (p,v), so that the line L becomes the new coordinate axis v = 0. Many features, that 
are symmetric with respect to L will then simplify; e.g. the distance to the bifurcation point along 
L is precisely u. 



a b\ ( /j, 
-b a)\v 



afj, + bv 
av — bji 



a —b 
b a 



(13) 



where a = sinaoi b — coscto- Figure [2j a) shows the relationship between the two coordinate 
systems. 



3.1 Fixed points and their local bifurcations 

The normal form ( ]12[ ), in terms of the modulus and phase z = re 1 *, is 

r = r(/i — ar 2 ) + cos (j), 

. „ 1 

cp = v — br sin <fi. 

r 



(14) 



The fixed points are given by a cubic equation in r 2 , and so we do not have convenient closed forms 
for the corresponding roots (the Tartaglia explicit solution is extremely involved). However, it is 
easy to obtain the locus where two of the fixed points coalesce. The parameter space is divided 
into two regions, region III with three fixed points, and the rest of parameter space with one fixed 
point, as seen in figure^lfo). The curve separating both regions is a saddle-node curve given by (see 







(3 + 3s 2 ,2V3s) 
(2 + 6s 2 ) 2 / 3 



s € (— oo, +oo), 



(15) 



in (u,v) coordinates (13). The saddle-noddle curve is divided into three different arcs SN + , SN_ 
and SNo by two codimension-two cusp bifurcation points, Cusp±. SN_ corresponds to values 
s e (—oo,l), SN to s € (— 1,+1) and SN + to s e (l,+oo). The cusp points Cusp± have values 
s = ±1. The curves SN + and SN_ are asymptotic to the line L, and region II is the pinning region 



G 




Figure 3: (a) Bifurcations of the fixed points corresponding to the normal form (12 1. SN± and SN 
are saddle-node bifurcation curves, H± are Hopf bifurcation curves, Cusp± are cusp bifurcation 
points and TB-t are Takens-Bogdanov bifurcation points. A perspective view of the corresponding 
codimension-thrcc bifurcation in terms of (/x, v 1 e) is shown in (b). (c) is a zoom of (a), showing that 
Cusp_ and TB_ are different. The parameter range in (a) is {v, jit) £ [—3,6] x [—1,6], for a = 45°. 



in this case. The two fixed points that merge on the saddle-node curve have phase space coordinates 
Z2 = ^e 1 ^ 2 , and the third fixed point is zq = roe 1 ^ , where 

r ° = liWJ ' r2 = l^H ' (16) 

and the phases are obtained from sin<p = r{y — br 2 ), cos</> = r(ar 2 — /i). 

The Hopf bifurcations of the fixed points can be obtained by imposing the conditions T = and 



D > 0, where T and D are the trace and determinant of the Jacobian of (12 1. These conditions 
result in two curves of Hopf bifurcations (see [B] for details): 

V a VI -s 2 J (17) 
(-l,-v/(l-6)/2), H+ : se (V(l + 6)/2,+l). 



H_ : s e 



For s — > ±1 both curves are asymptotic to the [i = axis [v — > ±oo), the Hopf curve for e = 0; 
the stable limit cycles born at these curves are termed C_ and C+ respectively. The limit cycles 
(rotating waves) in II± rotate in opposite directions, and III is the pinning region where the rotation 
stops and we have a stable fixed point. Solutions with u) = 0, that existed only along a single line in 
the absence of imperfections, now exist in a region of finite width. Figure [3^6) shows what happens 
when the e dependence is restored; what we have is that figure [3ja) just scales with e as indicated 



in (11 1, and the pinning region collapses onto the line L of the perfect case with 5*0(2) symmetry. 

The other ends of the H± curves are on the saddle-node curves of fixed points previously ob- 
tained, and at these points T = D = 0, so they are Takens-Bogdanov points TB±, as shown in 
figure [3j The TB_ and Cusp_ codimension-two bifurcation points are very close, as shown in the 
zoomed-in figure [3^6). In fact, depending on the angle ao, the Hopf curve H_ is tangent to, and 
ends at, either SN_ or SN . For a > 60°, H_ ends at SN , and for a = 60° Cusp_ and TB_ 
coincide, and H_ ends at the cusp point, a very degenerate case. 

From the Takens-Bogdanov points, dynamical systems theory says that two curves of homoclinic 
bifurcations emerge, resulting in global bifurcations around these points. Moreover, the stable limit 
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Figure 4: Schematic of the bifurcations of the normal form (12 1 
bifurcations 



There are seven curves of global 
(homoclinic collisions of a limit cycle with a saddle), CF (a cyclic- 



Hom±, Homo 

fold), SNIC±, and SNICo, and nine codimension-two points (black circles). The regions around the 
codimension-two points have been enhanced for clarity. 



cycles in regions do not exist in region III, so they must disappear in additional bifurcations. 
These additional global bifurcations have been explored numerically and using dynamical systems 
theory; these are summarized in figure [4j There are nine codimension-two points organizing the 
dynamics of the normal form (12). Apart from the cusp and Takens-Bogdanov points already 



found, Cusp± and TB±, there are five new points: Ba, CfHom and three different Snic- Homoclinic 
bifurcations, SnicHom± and SnicHomo. On H_, before crossing the SNo curve, the Hopf bifurcation 
becomes subcritical at the Bautin point Ba, and from this point a curve of cyclic-folds CF appears. 
This curve is the limit of the subcritical region where two periodic solutions, C_ and Co, exist 
and they merge on CF. Co is the unstable limit cycle born in the branch of H_ between Cusp_ 
and TB_, from now on termed the Hopf curve Ho. Inside the pinning region these two periodic 
solutions disappear when they collide with a saddle fixed point along the curves Homo (Co collision) 
and Hom_ (C_ collision), in the neighborhood of the Takens-Bogdanov point TB_, where the 
homoclinic curve Homo is born. Away from the TB_, for increasing values of [i and v, the curve 
Hom_ becomes tangent to and collides with the SN_ curve, in a SnicHom_ bifurcation. Very close 
to SnicHom_, SN_ closely followed by Hom_ become indistinguishable from the SNIC_ bifurcation, 
in the same sense as discussed in section §4.3[ 

The two curves Hom_ and Hom on approaching SN , result in a couple of codimension-two 
bifurcation points, SnicHom and CfHom. The arc of the curve SN between the two new points 
SnicHomo and CfHom, is a curve of saddle-node bifurcations taking place on the limit cycle Co, 
resulting in the SNICo bifurcation curve, as shown in figure|4j and in more detail in the numerically 
computed inset figure (5ja) . SnicHomo is exactly the same bifurcation as SnicHom±. 

The cyclic fold bifurcation curve CF intersects the SNICo bifurcation curve past the SnicHomo 
point, i.e. when the SN curve is a line of SNIC bifurcations, at the point CfHom. On the SNIC 
curve, one of the limit cycles born at CF undergoes a SNIC bifurcation. At the point CfHom, 
the SNIC bifurcation happens precisely when both limit cycles are born at CF: it is a saddle-node 
bifurcation of fixed points taking place on a saddle-node bifurcation of limit circles. After the 
CfHom point, the SNIC curve becomes an ordinary saddle-node bifurcation curve, and there is 
an additional homoclinic bifurcation curve emerging from this point CfHom, Hom_ in figure [4] 
The two limit cycles born at CF exist only on one side of the CF line, so when following a closed 
path around the CfHom point they must disappear. One of them undergoes a SNIC bifurcation on 
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1: ijl = 1.67732854 2: fi = 1.6776 3: fi = 1.6777322 




4: fi = 1.6784 5: fi = 1.678312 6: fi = 1.6782 7: (j, = 1.678028 8: fi = 1.6775 

Figure 5: (a) Zoom in parameter space around the codimension-two global bifurcation points 
SnicHom and CfHom. (6) Phase portraits around CfHom; plots 1 to 4 at v = 0.473, plots 5 to 8 at 
v = 0.475, for /j, as specified. 
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(a) 



(b) 



H Q B /SN_ 




HI SN + Po 




PF_ 



PF_ 



Bo -h - 



—ppr 



SN+ 



Figure 6: Steady bifurcations of fixed points corresponding to the normal form ( |20[ ): (a) regions in 
parameter space delimited by the fixed points and their steady bifurcations; (6) and (c) bifurcations 
along the paths A and B shown in (a), respectively. 



the SNIC curve and disappears. The other limit cycle collides with the saddle point born at the 
saddle-node curve SN (or SNICo) and disappears at the homoclinic collision Hom_. Numerically 
computed phase portraits illustrating these processes around the CfHom point are shown in figure [5] 
The other Takens-Bogdanov point TB + does not present any additional complications. The 
homoclinic curve Hom + emerging from it approaches and intersect the SN + curve in a SnicHom + 



the stable limit cycles in 



codimension two point, as shown in figure |4j For large values of /i 2 
regions II± disappear at SNIC± (saddle-node on an invariant circle) bifurcation curves. On these 
curves, a saddle-node bifurcation of fixed points takes place on top of the limit cycle, and the cycle 
disappears in an infinite-period bifurcation. What remains, and is observable, is the stable fixed 
point born at the saddle-node. 

The width of the pinning w(d, e) region away from the origin (large s) is easy to compute from 



(15) 



d = u~ (3s 2 /4) 1/3 , w = 2v ~ (16/V3s) 1/3 =*> w = 2/Vd. (18) 

Restoring the e dependence, we obtain w(d, e) = 2e/Vd. The pinning region becomes narrower 
away from the bifurcation point, and the width is proportional to e, the size of the imperfection. 



4 Symmetry breaking of SO (2) to Z2, with an ez term 



The term in (Joj) corresponds to breaking SO (2) symmetry in a way that leaves a system with 
Z2 symmetry, corresponding to invariance under a half turn. The normal form to be analyzed is 
(ffol with p = 1, q = and e = 1: 



z = z(/i + iv — c|z| 2 ) + z. 



(19) 



The new normal form ( 19 ) is still invariant to z — > —z, or equivalently, the half-turn <p — > (fy+n. This 



is all that remains of the SO (2) symmetry group, which is reduced to Z2, generated by the half-turn. 
In fact, the Z2 symmetry implies that P(z,z) in z = P(z,z) must be odd: P{— z,— z) = —P(z,z), 
which is ([5| for 8 = n, the half turn. Therefore, ( 19 ) is the unfolding corresponding to the symmetry 
breaking of SO{2) to Z 2 . This case has also been analyzed in Gambaudo] ( |1985| ; Wagener (2001 ); 
Broer et al.\ p008| ; [Saleh fc Wagener| ( [20To| . 

Writing the normal form ( 19 ) in terms of the modulus and phase of z = re 1 ^ gives 



f = r(/j, — ar 2 ) + r cos 20, 
(f> = v — br 2 — sin 2(f). 



(20) 



4.1 Fixed points and their bifurcations 



The normal form (191 
trivial solution Pq 



or (20), admits up to five fixed points. 



One is the origin r = 0, the 
The other fixed points come in two pairs of .^-symmetric points: one is the 
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pair P + = r + e 1< ^+ and P£ = — r + e Ic ^+, and the other pair is P_ = r^c 1 ^- and PI = — r_e"P-. 
Coefhcients r± and 0± are given by 

r\ = a f i + bv ± A, cf>± = (a ±a 1 )/2, (21) 
A 2 = 1 - (av- bp) 2 , e iai = av - bfi - iA. (22) 

The details of the computations are given in [C] There are three different regions in the (fi, v)- 
parameter plane: region III, where there exist five fixed points, Pq, P± and P±; region Io where 
three fixed poins exist, Pq, P + and PT; and the rest of the parameter space where only Pq exists. 
These three regions are separated by four curves along which steady bifurcations between the 
different fixed points take place, as shown in figure [6j Along the semicircle 

PF+ : fi 2 + v 2 = 1 and a^i + bv < 0, (23) 

the two symmetrically-related solutions P+ and P* are born in a pitchfork bifurcation of the trivial 
branch P . Along the semicircle 

PF_ : /r 2 + v 2 = 1 and api + bv > 0, (24) 

the two symmetrically-related solutions P_ and PI are born in a pitchfork bifurcation of the trivial 
branch Pq. Along the two half- lines 

SN+ : (i = (av - l)/b, SN_ : fj, = (av + l)/b, both with a/i + bv > 0, (25) 

a saddle-node bifurcation takes place. It is a double saddle-node, due to the Z 2 symmetry; we have 
one saddle-node involving P + and P_, and the ^-symmetric saddle-node between P* and P*. 
Bifurcation diagrams along the paths A and B in figure [6](a) are shown in parts (b) and (c) of the 
same figure. 

We can compare with the original problem with SO(2) symmetry, corresponding to e = 0. In 
order to do that, the e dependence will be restored in this paragraph. The single line L (jj, = v tan ag) 
where uj = and nontrivial fixed points exist in the perfect problem, becomes a region of width 2e 
in the imperfect problem, where up to four fixed points exist, in addition to the base state Pq; they 
are the remnants of the circle of fixed points in the original problem. Solutions with uj = 0, that 
existed only along a single line in the absence of imperfections, now exist in a region bounded by the 
semicircle P + and the half-lines av — bfi = ±e; this region will be termed the pinning region. It bears 
some relationship with the frequency-locking regions appearing in Neimark-Sacker bifurcations, in 
the sense that here we also have frequency locking, but with lu = 0. The width of the pinning region 
is proportional to e, a measure of the breaking of 50(2) symmetry due to imperfections. 

In the absence of imperfections (e = 0) the Pq branch looses stability to a Hopf bifurcation 
along the curve fj, = 0. Let us analyze the stability of Pq in the imperfect problem. Using Cartesian 



coordinates z = x + iy in ( 19 ) we obtain 



iM"" 1 fl -i)(;)- (i2+ » 3 »fe;a::)- <*> 



The eigenvalues of Pq are the eigenvalues of the linear part of ( 26 1 , X± — /i ± \/l — v 2 . There is 
a Hopf bifurcation (^s\± ^ 0) when /i = and \v\ > 1, i.e. on the line /i = outside region II; 
the Hopf frequency is uj — sign(zy) v / ^ 2 — 1- The sign of ui is the same as the sign of v, from the 4> 



equation in |20j). Therefore, we have a Hopf bifurcation with positive frequency along H + (/i = 
and v > 1) and a Hopf bifurcation with negative frequency along H_ (/x = and v < 1). The 
bifurcated periodic solutions are stable limit cycles C+ and C_, respectively. 

The Hopf bifurcations of the P± and Pjjl points can be studied analogously. The eigenvalues of 



the Jacobian of the right-hand side of ( 26 1 at a fixed point characterize the different bifurcations 
that the fixed point can undergo. Let T and D be the trace and determinant of J. The eigenvalues 
are given by 

A 2 -TA + D = => A = -(T± V / Q), Q = T 2 - AD. (27) 
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Figure 7: Local bifurcations of fixed points in the symmetry breaking of SO(2) to Z2 case. 
Codimension-one bifurcation curves: Hopf curves H± and Ho, pitchforks PF±, and saddle-nodes 
SN±. Codimension-two bifurcation points: degenerate pitchforks dPF±, Takens-Bogdanov TB and 
TB±. L is the zero-frequency curve in the SO(2) symmetric case. 



A Hopf bifurcation takes place for T — and Q < 0. The equation T — at the four points P± 
results in the ellipse /1 2 — 4afe/i^ + 4a 2 v 2 — 4a 2 (see [c] for details). This ellipse is tangent to SN_ 
at the point (p,, v) = (2b, (b 2 — a 2 )/a). The condition Q < is only satisfied by P + and PI on the 
elliptic arc H from (fj,,v) — (0,-1) to (26, (b 2 — a 2 ) /a): 

fj, = 2abv + 2ay/l - a 2 v 2 , v € [-1, (b 2 - a 2 )/a\. (28) 

The elliptic arc Ho is shown in figure [7j Along this arc a pair of unstable limit cycles Co and Cg 
are born around the fixed points P+ and P£, respectively. 



4.2 Codimension-two points 

The local codimension-one bifurcations of the fixed points are now completely characterized. There 
are two curves of saddle-node bifurcations, two curves of pitchfork bifurcations and three Hopf bi- 
furcation curves. These curves meet at five codimension-two points. The analysis of the eigenvalues 
at these points, and of the symmetry of the bifurcating points (Pq, P + and PI), characterizes these 
points as two degenerate pitchforks dPF±, two Takens-Bogdanov bifurcations with Z2 symmetry 
TB±, and a double Takens-Bogdanov bifurcation TB, as shown in figure [7| 

The degenerate pitchforks dPF± correspond to the transition between supercritical and subcrit- 
ical pitchfork bifurcations. At these points, a saddle-node curve is born, and only fixed points are 
involved in the neighboring dynamics. The only difference between dPF + and dPF_ is the stability 
of the base state Po; it is stable outside the circle fi 2 + v 2 = 1 at dPF + and unstable at dPF_. 
Schematics of the bifurcations along a one-dimensional path in parameter space around the dPF + 
and dPF_ points are illustrated in figure [8] The main difference, apart from the different stability 
properties of Pq, P + and P^_, is the existence of the limit cycle C_ surrounding the three fixed 
points in case (b), dPF_. 



The Takens-Bogdanov bifurcation with Z2 symmetry has two different scenarios ( Chow et al. 



19941, and they differ in whether one or two Hopf curves emerge from the bifurcation point. In our 
problem, bifurcation point TB + has a single Hopf curve, H + , while the TB_ point has two Hopf 
curves, H_ and Ho, emerging from the bifurcation point. The scenario TB + is depicted in figurc[9j 
showing the bifurcation diagram as well as the bifurcations along a closed one-dimensional path 
around the codimension-two point. We have also included the P+ and PI solutions that merge with 
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B c » 1 { 1 oB 

PF+ PF_ SN_ 



Figure 8: Schematics of the degenerate pitchfork bifurcations (a) dPF + and (b) dPF_. On the 
left, bifurcation curves emanating from dPF± in parameter space are shown, along with a closed 
one-dimensional path (dashed); shown on the right are schematics of the bifurcations along the 
closed path, starting and ending at A (dPF + ) and B (dPF_). The fixed point curves are labeled 
with the signs of their eigenvalues. C_ is the periodic solution born at the curve H_ . 
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a: fj, = 0.05, v = 0.9 b: fj, = 0.1, i/ = 1.1 




Figure 9: Takens-Bogdanov bifurcation with Z 2 symmetry TB + . The top left shows bifurcation 
curves emanating from TB + in parameter space, with a closed one-dimensional path. The top right 
shows a schematic of the bifurcations along the closed path, starting and ending at A. The fixed 
point curves are labeled with the signs of their eigenvalues. C+ is the periodic solution born at the 
curve H + . The region inside the dashed rectangle on the right contains the states locally connected 
with the bifurcation TB + . The bottom panels show four numerically computed phase portraits, at 
points labeled a, b, c and d, for the specified parameter values. 
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Figure 10: Takens-Bogdanov bifurcation with Z2 symmetry, TB_. The top left shows bifurcation 
curves emanating from TB_ in parameter space, with a closed one-dimensional path. The top right 
shows a schematic of the bifurcations along the closed path, starting and ending at A. The fixed 
point curves are labeled with the signs of their eigenvalues. C_ is the periodic solution born at 
the curve H ; C + and Cq + are the unstable cycles born simultaneously at the Hopf bifurcation 
Ho around the fixed points P+ and PT; and Co is the cycle around both fixed points that remains 
after the gluing bifurcation. The bottom panels show three numerically computed phase portraits, 
at points labeled a, b and c, for [i = 0.5 and v as indicated, illustrating the gluing and cyclic fold 
bifurcations. 



the P_ and PI fixed points along the saddle-node bifurcation curve SN + , although they are not 
locally connected to the codimension-two point, in order to show all the fixed points in the phase 
space of (19). A curve of global bifurcations, a heteroclinic cycle Het + connecting P_ and PI, is 
born at TB + . The heteroclinic cycle is formed when the limit cycle C+ simultaneously collides with 
the saddles P_ and P* . 



The scenario TB_, is depicted in figure 10 showing the bifurcation diagram and also the bi- 
furcations along a closed one-dimensional path around the codimension-two point. Two curves of 
global bifurcations are born at TB_. One corresponds to a gluing bifurcation Glu, when the two 
unstable limit cycles C + and CJ + , born at the Hopf bifurcation Ho around P+ and PT , simultane- 
ously collide with the saddle Po; after the collision a large cycle Co results, surrounding the three 
fixed points P , P + and P*. The second global bifurcation curve corresponds to a saddle-node of 
cycles, where Co and C_ collide and disappear. 

A generic Takens-Bogdanov bifurcation (without symmetry) takes place at the TB point on 
the SN_ curve. At the same point in parameter space, but separate in phase space, two Z2 sym- 
metrically related Takens-Bogdanov bifurcations take place, with P + and PT being the bifurcating 
states. A schematic of the bifurcations along a one-dimensional path in parameter space around 
the TB point is shown in figure [TT] Apart from the states locally connected to both TB bifurcation, 
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there also exist the base state P Q and the limit cycle C_ 



4.3 Global bifurcations 

In the analysis of the local bifurcations of fixed points we have found three curves of global bifur- 
cations, a gluing curve Glu and a saddle-node of cycles CF emerging from TB_, a heteroclinic loop 
born at TB + (Het+), and a homoclinic loop emerging from TB (Horn). One wonders about the 
fate of these global bifurcation curves, and about possible additional global bifurcations. Numerical 
simulations of the solutions of the normal form ODE system (|T9| , or equivalently ( 26 ) , together 



with dynamical systems theory considerations have been used to answer these questions, and a 



schematic of all local and global bifurcation curves is shown in figure 12 

The gluing bifurcation Glu born at TB_ and the two homoclinic loops emerging from the two 
Takens-Bogdanov bifurcations TB (bifurcations of the symmetric fixed points P+ and PI) meet at 
the point PfGl on the circle PF_, where the base state Pq undergoes a pitchfork bifurcation (see 
figure [llfo). At that point, the two homoclinic loops of the gluing bifurcation, both homoclinic 
at the same point on the stable Po branch, split when the two fixed points P_ and P* bifurcate 
from P (see figure 13 The two homoclinic loops are then attached to the bifurcated points and 
separate along the curve Horn. The large unstable limit cycle Co, after the pitchfork bifurcation 
PF_ , collides simultaneously with both P_ and PI , forming a heteroclinic loop along the curve Hcto 
(see figure 13: and d2). Both curves Horn and Heto are born at PfGl and separate, leaving a region 
in between where none of the cycles C +, Cp + and Co exist. The unstable periodic solution Co 
merges with the stable periodic solution C_ that was born in H_ and existed in region I, resulting 
in a cyclic-fold bifurcation of periodic solutions CF. Phase portraits around PfGl are shown in 



figure 13 [d). 

The curve Horn born at PfGl ends at the double Takens-Bogdanov point TB. Locally, around 
both Takens-Bogdanov bifurcations at TB, after crossing the homoclinic curve the limit cycles C + 
and Cp + disappear, and no cycles remain. The formation of a large cycle Co at Heto surrounding 
both fixed points P+ and P+ is a global bifurcation involving simultaneously both P_ and Pt 
unstable points. It is the re-injection induced by the presence of the Z2 symmetry that is responsible 
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bifurcation curves Heto an d CF become very close when leaving the PfGl neighborhood, and merge 
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Figure 12: Global bifurcations in the symmetry breaking of 50(2) to case. Codimension-one 
bifurcation curves: Gluing Glu, cyclic-fold CF, homoclinic collision Horn and heteroclinic loops 
Heto, Het±. Codimension-two bifurcation points: pitchfork-gluing bifurcation PfGl, cyclic-fold 
heteroclinic bifurcation CfHet and two SNIC-heteroclinic bifurcations SnicHet±. 



at some point in a CfHom (Cyclic-Fold-Heteroclinic collision) codimension-two global bifurcation. 
After CfHom, the stable limit cycle C_, instead of undergoing a cyclic-fold bifurcation, directly 
collides with the saddle points P_ and P*_ along the Het_ bifurcation curve (see figure 12 1. In fact 
Heto and Het_ are collisions of a limit cycle with P_ and P*, but the limit cycle is on a different 
branch of the saddle-node of cycles CF on each side of CfHom. The two limit cycles born at CF 
are extremelly close together in the neighborhood of CfHom, and it is impossible to see them in a 
phase portrait, except with a very large zoom around P_ or P*. 

When increasing fi 2 + v 2 , the heteroclinic loop Het_ born at CfHet intersects the SN_ curve at 
a codimension-two global bifurcation point SnicHct'. When they intersect, the saddle-node appears 
precisely on the limit cycle, resulting in a SNIC bifurcation (a saddle-node on an invariant circle 
bifurcation). At the SnicHet_ point, the saddle-node and homoclinic bifurcation curves become 
tangents, and the saddle- node curve becomes a SNIC bifurcation curve for larger values of /i 2 + v 2 . 
Figure 14 shows numerically computed phase portraits for olq — 45°, below and above the SnicHet_, 
located at v « 0.1290, fi ~ 1.505. In the first case (figure 14 z) the saddle-node bifurcation SN_ 
takes place in the interior of C_, while in the second case (figure 14)) SN_ happens precisely on 
top of C_, resulting in a SNIC bifurcation. 

There remains a global bifurcation curve to be analyzed, the heteroclinic loop Het born at 
TB+. As shown in figure [l2j the curve Het intersects the SN + curve tangentially at a codimension- 
two global bifurcation point SnicHet. Beyond this point, the SN + curve becomes a line of SNIC 
bifurcations, where the double saddle- node bifurcations appear on the stable limit cycle C+, which 
disappears on entering the pinning region III, exactly in the same way as has been discussed for 
the SnicHet' bifurcation. 

We can estimate the width of the pinning region as a function of the magnitude of the imperfec- 
tion e and the distance d to the bifurcation point, w(d, e). The distance d will be measured along 
the line L, and the width w(d) will be the width of the pinning region measured transversally to L 
at a distance d from the origin. It is convenient to use the coordinates (u, v), so that the parameter 
u along L is precisely the distance d. If the pinning region is delimited by a curve of equation 
v = ±h(u), then w — 2h(u) — 2h(d). With an imperfection of the form ez, the case analyzed in 
this section, the pinning region is of constant width w = 2. By restoring the dependence on e, we 
obtain a width of value w(d, e) — 2e, independent of the distance to the bifurcation point; the width 
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Figure 13: Pitchfork-gluing bifurcation PfGl. Top left, bifurcation curves emanating from PfGl in 
parameter space; top right, schematics of the bifurcations along the line Glu-Hom. The bottom row 
shows four numerically computed phase portraits, at points labeled 1, 2, 3 and 4, for the specified 
parameter values. 
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(a) = 1.6, v = 0.19 (6) /z = 1.6, v = 0.213 




Figure 14: Numerically computed phase portraits in the ez case, for a = 45° and fi and v as 
indicated; cases (a) and (6) are below the SnicHet_ point and cases (c) and (d) are above the 
SnicHct_ point. 
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Figure 15: Bifurcation curves corresponding to the normal forms with quadratic terms (291 in the 
ez 2 case, (a) ao > 30°, (b) ao < 30°. Ho is tangent to the parabola at the Takens bogdanov points 
TB±, and H , Hom almost coincide with SN_; in the figure the distances have been exagerated 
for clarity. 

of the pinning region is proportional to e. 

5 Symmetry breaking of 50(2) with quadratics terms 

5.1 The ez 2 case 

The normal form to be analyzed in this case is 

i = z(fi + iv - c\z\ 2 ) + z 2 , (29) 
or in terms of the modulus and phase z = re 1 ^, 

r = r(ri — ar 2 ) + r 2 cos3<i, 

• o (30) 

<j) = v — br~ — r sin 

The fixed points are the origin Pq (r = 0) and the solutions of the same biquadratic equation as in 
the two previous cases. However, there is an important difference: due to the factor 3 inside the 



trigonometric functions in ( 30 1, the P_j_ points come in triplets (i = 1,2,3), each triplet has the same 



radius r but their phases differ by 120°. This is a consequence of the invariance of the governing 



equation (29) to the Z 3 symmetry group generated by rotations of 120° around the origin. This 



invariance was not present in the two previous cases. The bifurcation curves of the fixed points 



(excluding the Hopf bifurcations of P_j_ ) are still given by figure 17 a), but now a triplet of symmetric 
saddle-node bifurcations take place simultaneously on the SN + and SN_ curves. 

It can be seen (details inp3| that P % _ are saddles in the whole of region III, while Pj_ are stable, 



except for small angles ao < 7r/6 in a narrow region close to SN_ where they are unstable. For 



ao > it/ 6, the bifurcation diagram is exactly the same as in the zz case (figure 16 1), except that 
the homoclinic curves are now heteroclinic cycles between the triplets of saddles PI; this case is 
illustrated in figure [l5[a) . For ao < 7r/6, two Takens-Bogdanov bifurcation points appear at the 



tangency points between the SN_ curve and the arc H of the ellipse (b/i — 2av) + (afi — 1) = 1 



as shown in figure 15 b). The arc Ho is a Hopf bifurcation curve of P]_: three unstable limit 
cycles Cg are born when the triplet Pj. becomes stable. These unstable limit cycles disappear 
upon colliding with the saddles PI on a curve of homoclinic collisions, Homo, that ends at the 
two Takens-Bogdanov points TB + and TB_. This situation is very similar to what happens in 
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the ez case analyzed in section Sj4j where a Hopf curve H appeared close to SN_ joining two 
Takens-Bogdanov points. In both cases, the SO(2) symmetry is not completely broken, but a Z m 
symmetry remains. For «o — tt/6, the ellipse Ho becomes tangent to SN_ and the two Takens- 
Bogdanov points coalesce, disappearing for ao > tt/6. 

Finally, the stable limit cycles C_ and C+, born at the Hopf bifurcations H_ and H + , upon 
entering region III collide simultaneously with the saddles Pi, i = 1, 2 and 3, and disappear along 
two heteroclinic bifurcation curves Het and Hct' for small values of [i. These curves collide with 
the parabola at the codimension-two bifurcation points SnicHct and SnicHct'. For larger values of 
[i the limit cycles C_ and C+ undergo SNIC bifurcations on the parabola. The two curves Hom± 
emerge from the origin, as in the previous cases. 

In the three quadratic cases, the pinning region is delimited by the same parabola u = v 2 — 1/4, 
using the (u, v) coordinates introduced in ( |l3| ) (see also figure[2ja)). The width of the pinning region 
is easy to compute, and is given by w — 2v — 2^/u + 1/4 ~ 2\fd. By using (111, the dependence 



on e is restored, resulting in w(d,e) = 2eVd. The width of the pinning region increases with the 
distance d to the bifurcation point, and it is proportional to the amplitude of the imperfection e. 

5.2 The ezz case 

The normal form to be analyzed in this case is 

i = z(p + iv — c\z\ 2 ) + zz, (31) 
or in terms of the modulus and phase z — re 1 ^, 

r = 7'(/i — ar 2 ) + r 2 cos (f), 



= v — br — 



., . (32) 



r sin ( 



There are three fixed points: the origin P ( r = 0) and the two solutions P± of the biquadratic 
equation r 4 — 2(a/i + bv + l/2)r 2 + /i 2 + v 2 = 0, given by 



r± = a[i + bv + - ± y a^i + bv + - — (av — b[i) 2 . (33) 

These solutions are born at the parabola a/i + bv + 1/4 = (av — fr/i) 2 , and exist only in its interior, 
which is the pinning region III in figure [l6j The parabola is a curve of saddle-node bifurcations. It 
can be seen (details in[D]) that P+ is stable while P_ is a saddle in the whole of region III, so there 
are no additional bifurcations of fixed points in the zz case, in contrast to the z 2 case. 

As the perturbation is of second order, the Jacobian at Pq is the same as in the unperturbed case, 
and the Hopf bifurcations of Pq take place along the horizontal axis = 0. As in the unperturbed 
case, the Hopf frequency is negative for v < (H_), it is zero at the origin, and becomes positive for 
v > (H + ). The Hopf curves H_ and H + extend in this case up to the origin, in contrast with the 
zero and first order cases examined previously, where the Hopf curves ended in Takens-Bogdanov 
bifurcations without reaching the origin. 

The stable limit cycles C_ and C+, born at the Hopf bifurcations H_ and H + , upon entering 
region III collide with the saddle P+ and disappear along two homoclinic bifurcation curves Hom± 
for small values of /i. For larger values of /i, the curves Hom± collide with the parabola at the 
codimension-two bifurcation points SnicHom-t, and for larger values of /i, the saddle-node bifurca- 
tions take place on the parabola and the limit cycles C_ and C + undergo SNIC bifurcations, as in 
the previous e and ez cases. Figure [l6^ a) summarizes all the local and global bifurcation curves in 
the zz case, and shows numerically computed phase portraits around the SnicHom_ point. Fig- 
ure 16 6) shows the SN + and Hom_ bifurcations before SnicHom_ (at /i = 0.03), and figure 16 'c) 



illustrates the SNIC_ bifurcation after SnicHom_ (at fi = 0.033). 
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v = -0.306 v = -0.30442 v = -0.303 



(c) 




v = -0.307 v = -0.30552 v = -0.304 



Figure 16: (a) Schematic of bifurcation curves corresponding to the normal form with quadratic 
terms in the ezz case. Phase portraits (b) crossing the Hom_ curve at /i = 0.03, for v values 
as specified, and (c) crossing the SNIC_ curve at fi = 0.033. Thick grey lines correspond to the 
periodic orbit and the homoclinic and heteroclinic loops. 
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(a) 



SNIC_ 




SNIC_ 



SnicHom^ 




Figure 17: (a) Bifurcation curves corresponding to the normal form with quadratic terms in the 
ez 2 case, (b) Phase portrait on the curve H showing the degenerate Hopf bifurcation and the 
associated homoclinic loop. 



5.3 The ez 2 case 



The normal form to be analyzed in this case is ( 10 ) with p = q = 2 and e = 1: 

i = + iv — c|z| 2 ) + z 2 . 



The normal form ( 34 ) , in terms of the modulus and phase z = re 1 ^ , reads 



r = r(/i — ar 2 ) + r 2 cos (j), 
6 = v — br 2 + r sin d>. 



(34) 



(35) 



There are three fixed points: the origin Po (r = 0) and the two solutions P± of the biq uadr atic 
equation r 4 — 2(a/i + bv + l/2)r 2 + /.i 2 + v 2 = 0, which are the same as in the ezz case (( 5.2 1. In 



fact, the fixed points in these two cases have the same modulus r and their phases have opposite 
sign; changing <f> — > —0 in (32) results in (35). Therefore, the bifurcation curves of the fixed points 



(excluding Hopf bifurcations of P±) in this case are also given by figure [l6[a). 

It can be seen (details in|D| that P+ is a saddle in the whole of region III, while P_ is stable 
for fi > 0, unstable for fi < 0, and undergoes a Hopf bifurcation Ho along the segment of fi = 
delimited by the parabola of saddle-node bifurcations. The points TB-t where Ho meets the parabola 
are Takens-Bogdanov codimension-two bifurcations. Figure [17] summarize all the local and global 
bifurcation curves just described. 

In the present case, the two Takens-Bogdanov bifurcations and the Hopf bifurcations along Ho 
are degenerate, as shown in [D] Detailed analysis and numerical simulations show that the Hopf 
and homoclinic bifurcation curves emerging from the Takens-Bogdanov point are both coincident 
with the Ho curve previously mentioned. Moreover, the interior of the homoclinic loop is filled with 
periodic orbits, and no limit cycle exists on either side of H . This situation is illustrated in the 



phase portrait in figure 17 ^6) . This highly degenerate situation will be broken by the presence of 



additional terms in the normal form, and of the continuous family of periodic orbits, only a few will 
remain. Dumortier et al. (1987), who have analyzed in detail the unfolding of such a degenerate 



case, find that at most two of the periodic orbits survive. 
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(a) m = 4 



(6) m > 4 




Finally, and exactly in the same way as in the zz case examined in the previous subsection, the 
stable limit cycles C_ and C+ born respectively at the Hopf bifurcations H_ and H + and existing 
in regions I_ and I + (figure 17 z), on entering region III collide with the saddle P + and disappear 
along two homoclinic bifurcation curves Homj- for small values of fi. These curves collide with the 
parabola at the codimension-two bifurcation points SnHom-t and for larger values of /z, the limit 
cycles C_ and C + undergo SNIC bifurcations on the parabola. The two curves Hom-t emerge from 
the origin, as in the previous case. 



6 Symmetry breaking 50(2) — > Z m , m > 4 

For completeness, and also for intrinsic interest, we will explore the breaking of the 50(2) symmetry 
to Z m , so that the imperfections added to the normal form ^ preserve the Z m subgroup of SO (2) 
generated by rotations of 2w/m (also called C m ). The lowest order monomial in (z,z) not of the 
form z\z\ 2p and equivariant under Z m , is z m_1 , resulting in the normal form 



i = z(fi + if — c|z| 2 ) + ez 



m — l 



(36) 



In terms of the modulus and phase of the complex amplitude z = re 1 ^, the normal form becomes 

r = r(n - ar 2 ) + er TO_1 



= v - br 2 - er m - 2 



cosmq 
sin mch. 



(37) 



= 2 and m = 3 have already been examined in Qand ^5.1 respectively. When m > 4, 
_1 is smaller than the remaining terms in (36 1, so the effect of the symmetry breaking 



The cases m 
the term ez" 

is going to be small compared with the other cases analyzed in this paper. The fixed point solutions 
of (37), apart from the trivial solution Pq (r = 0), are very close the zero-frequency line L in the 



points of (37) satisfy 



perfect system. Using the coordinates (it, v) along and orthogonal to L (13), the nontrivial fixed 

(r 2 



uf 



e 2 r 2m- 



(38) 



On L, r 2 — u; close to L, the fixed points P± are given by r 2 ~ u ± \J e 2 u m 2 — v 2 . The pinning 



at dominant order in e, is v = ±eu^ 



-2)/2 



region 

m — 4 and a horn for m > 4, as illustrated in figure 18 



This gives a wedge-shaped region around L for 
On the boundaries of the pinning region, 
the fixed points merge in saddle-node bifurcations that take place on the limit cycles C±. These 
are curves of SNIC bifurcations, exactly the same phenomena that is observed in Neimark-Sacker 



bifurcations (Arrowsmith & Place 1990), and that we have encountered also in the previous cases 



analyzed in the present study. Due to the symmetry Z m , from (37) we see that at the boundaries 



of the wedge, m simultaneous saddle-node bifurcations take place. Figure 19 shows how the fixed 
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(a) 



(b) 




Figure 19: Crossing the horn for m — 5. (a) Fixed point solutions of (37) at the five points in 
figure [l8^£>); grey points in 2 and 4 are the saddle-node points, that split in a stable point (black) 
and saddle (white), (b) Phase portraits corresponding to the five cases in (a). 



points appear and disappear in saddle-node bifurcations on the limit cycle C± when crossing the 
horn for the m — 5 case, at the five points in parameter space indicated in figure 18 6). The 



nontrivial fixed points, from (37) and r ~ [i/a ~ v/b, satisfy 



r 2 = fj,/a + e(/i/a) (m 1)/2 costo</), 
r 2 = v/b- e(v/b) {m - 1)/2 shim(/), 



(39) 



and the intersection of these circles modulated by the sin mip and cos m</> terms is illustrated in fig- 
ure|19[a). Phase portraits corresponding to the five points in parameter space are also schematically 



shown in figure 19 b) 



The width of the pinning region in the symmetry breaking SO(2) — > Z m case is obtained from 
the shape of the horn region and is given by w(d, e) = 2v = 2eS m ~ 2 ^ 2 . Again, as in all the preceding 
cases, the width of the pinning region is proportional to the amplitude of the imperfection e. 



7 Common features in the different ways to break SO(2) 
symmetry 

Here we summarize the features that are common to the different perturbations analyzed in the 
previous sections. The most important feature is that the curve of zero frequency splits into two 
curves with a region of zero- frequency solutions appearing in between (the so-called pinning region). 
Of the infinite number of steady solutions that exist along the zero-frequency curve in the perfect 
system with 50(2) symmetry, only a small finite number remain. These steady solutions correspond 
to the pinned solutions observed in experiments and in numerical simulations, like the ones to be 
described in SjSJ The number of remaining steady solutions depends on the details of the symmetry- 
breaking imperfections, but when 5*0(2) is completely broken and no discrete symmetries remain, 



there are three steady solutions in the pinning region III (see figure 20 1) . One corresponds to 
the base state, now unstable with eigenvalues (+, +). The other two are born on the SNIC curves 
delimiting region III away from the origin. Of these two solutions, one is stable (the only observable 



state in region III) and the other is a saddle (see figure 20) and c). There are also the two Hopf 
bifurcation curves H_ and H + . The regions where the Hopf bifurcations meet the infinite-period 
bifurcations cannot be described in general, and as has been shown in the examples in the previous 
sections, will depend on the specifics of how the 50(2) symmetry is broken, i.e. on the specifics of 
the imperfections present in the problem considered. These regions contain complex bifurcational 
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Figure 20: Imperfect Hopf under general perturbations: (a) regions in parameter space; (b) and (c), 
bifurcation diagrams along the two one-dimensional paths, (1) and (2), respectively. The signs (++, 
. . . ) indicate the sign of the real part of the two eigenvalues of the solution branch considered. 
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Figure 21: Square root and logarithmic fits to the periods of C_ approaching the Hom_ curve 
at v = 0.6: (a) fitted curves in the complete range fi 6 [2.0135,2.2], and (b) fits in the range 
(i £ [2.0135, 2.027] for the log fit and fx € [2.018, 2.2] for the square root fit. 



processes, and are represented as grey disks in figure 20 ^o) . The stable limit cycle existing outside 



III, in regions I±, undergoes a SNIC bifurcation and disappears upon entering region III (see 
figure [20p) . When the SNIC bifurcation curves approach the Hopf bifurcation curves (i.e. enter the 
grey disks regions), the saddle-node bifurcations do not occur on the stable limit cycle but very 
close, and the limit cycle disappears in a saddle-loop homoclinic collision that occurs very close to 
the saddle-node bifurcations. These homoclinic collisions behave like a SNIC bifurcation, except in 
a very narrow region in parameter space around the saddle-node curves, as has been discussed in 

In all cases considered, the width of the pinning region scales linearly with the strength of 
the symmetry breaking e. In all cases, we have found w(d,c) = 2eS p ~ 1 ^ 2 , where p is the order 
of the symmetry breaking considered. For lower order terms, the width decreases (e case, order 
zero) or remains constant (ez case, order one) with increasing distance from the bifurcation point. 
For quadratic and higher order terms, the width increases with the distance. When arbitrary 
perturbations are included, we expect a behavior of the form w(d,e) = ef(d), where the function 
/ will depend on the details of the symmetry-breaking terms involved. The size of the regions 
containing complex bifurcational processes (the grey disks in figure 20 1) is of order e or smaller, as 
we have seen in all cases considered. Therefore these regions are comparable in size or smaller than 
the width of the pinning region. 



7.1 SNIC versus homoclinic collision: scaling of the period 

Another feature we have found in the five scenarios discussed above is that the SNIC± bifurcations, 
where the stable limit cycle disappears on entering the pinning region away from the origin (p = v = 
0), become homoclinic (or hctcroclinic if the SO (2) symmetry is not completely broken) collisions 
between the stable limit cycle and the saddle point that is born in the saddle-node bifurcations 
SN-t that now take place not on the limit cycle but very close to it. There is a codimension-two 
global bifurcation, termed SnicHom in the bifurcation diagrams discussed above, where the SNIC 
curve, the saddle-node curve and the homoclinic collision curve meet. 

The scaling laws of the periods when approaching a homoclinic or a SNIC bifurcation are 
different, having logarithmic or square root profiles: 

^Het = T m " 1" 0(1), T SNIC = — J=+0(1), (40) 

A fi - ^ VA 1 - 
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where A is the positive eigenvalue of the saddle, and k a constant. So the question is what happens 
with the scaling of the period close to the SnicHom bifurcation? 

We have analized in detail the ez case, being the other cases very similar. We refer to the 



figures 12 and 14 From a practical point of view, close to but before the SnicHet_ point, the saddle- 
node bifurcation SN_ is very closely followed by the heteroclinic collision of the limit cycle C_ 
with the saddles P_ and P*, and it becomes almost indistinguishable from the SNICL bifurcation. 
We have numerically computed the period of C_ at v = 0.6 for decreasing /i values approaching 
SN_ in the range fj, € [2.0135,2.2], for ao = 45°. Figure [2l|a) shows both fits using the values 
of the period over the whole computed range. The log fit overestimates the period while the 
square-root fit underestimates it, and this underestimate gets larger as the heteroclinic collision 



is approached. Figure 21 b) again shows both fits, but now using values close to the collision for 
the log fit and values far away from the collision for the square-root fit. Both fits are now very 
good approximations of the period in their corresponding intervals, and together cover all the values 
numerically computed. When the interval between the SN_ bifurcation and the heteroclinic collision 



(in figure 21 ^sn = 2.01420 and fi c — 2.01336, respectively) is very small, it cannot be resolved 
experimentally (or even numerically in an extended systems with millions of degrees of freedom, as 
is the case in fluid dynamics governed by the three-dimensional Navier-Stokes equations). In such 
a situation the square-root fit looks good enough, because away from the SN_ point, the dynamical 
system just feels the ghost of the about to be formed saddle-node pair and does not distinguish 
between whether the saddle-node appears on the limit cycle or very close to it. However, if we 
are able to resolve the very narrow parameter range between the saddle-node formation and the 
subsequent collision with the saddle, then the log fit matches the period in this narrow interval 
much better. 

Due to the presence of two very close bifurcations (Hef and SN_), the scaling laws become cross- 
contaminated, and from a practical point of view the only way to distinguish between a SNIC and 
a Homoclinic collision is by computing or measuring periods extremelly close to the infinite-period 
bifurcation point. We can also see this from the log fit equation in ( |40"| ); when both bifurcations are 
very close, A, the positive eigenvalue of the saddle, goes to zero (it is exactly zero at the saddle-node 
point), so the log fit becomes useless, except when the periods are very large. 

When the interval between the SN_ bifurcation and the homoclinic collision is very small, it 
cannot be resolved experimentally (or even numerically in an extended systems with millions of 
degrees of freedom, as is the case in fluid dynamics governed by the three-dimensional Navier- 
Stokes equations). In such a situation, the square-root fit appears to be good enough, because 
away from the SN_ point, the dynamical system just feels the ghost of the about-to-be-formed 
saddle-node pair and does not distinguish between whether the saddle-node appears on the limit 
cycle or just very close to it. However, if we are able to resolve the very narrow parameter range 
between the saddle-node formation and the subsequent collision with the saddle, then the log fit 
matches the period in this narrow interval much better. 

7.2 Codimension-two bifurcations of limit cycles 

The bifurcations that a limit cycle can undergo have been an active subject of research since 
dynamical systems theory was born. Even in the case of isolated codimension-one bifurcations, a 



complete classification was not completed until fifteen years ago (Turaev & Shilnikov 19951, when 
the blue-sky catastrophe was found. The seven possible bifurcations are: the Hopf bifurcation, 
where a limit cycle shrinks to a fixed point, and the length of the limit cycle reduces to zero. Three 
bifurcations where both the length and period of the limit cycle remain finite: the saddle node of 
cycles (or cyclic fold), the period doubling and the Neimark-Sacker bifurcations. Two bifurcations 
where the length remains finite but the period goes to infinity: the collision of the limit cycle with 
an external saddle forming a homoclinic loop, and the appearance of a saddle-node of fixed points 
on the limit cycle (the SNIC bifurcation). Finally, there is the blue-sky bifurcation where both 
the length and period go to infinity, corresponding to the appearance of a saddle-node of limit 
cycles transverse to the given limit cycle. The seven bifurcations are described in many books 



on dynamical systems, e.g. Shil'nikov et al. (2001); Kuznetsov (2004); they are also described on 
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the web page http://www.scholarpedia.org/article/Blue-sky_catastroplie maintained by A 



Shil'nikov and D. Turaev. 

Of the seven bifurcations, only four (Hopf, cyclic fold, homoclinic collision and SNIC) are 
possible in planar systems, as is the case in the present study, and we have found the four of them 
in the different scenarios explored. We have also found a number of codimension-two bifurcations 
of limit cycles. For these bifurcations a complete classification is still lacking, and it is interesting 
to list them because some of the bifurcations obtained are not very common. The codimension-two 
bifurcations of limit cycles associated to codimension-two bifurcations of fixed points can be found 
in many dynamical systems texbooks, and include Takcns-Bogdanov bifurcations (present in almost 
all cases considered here) and the Bautin bifurcation (in the e case). Codimension-two bifurcations 
of limit cycles associated only to global bifurcations are not so common. We have obtained five of 
them, that we briefly summarize here. 

PfGl A gluing bifurcation with the saddle point undergoing a pitchfork bifurcation; it may happen 
in systems with Z2 symmetry; see 

CfHom A cyclic-fold and a homoclinic collision occurring simultaneously; see ^4] 

CfHet A cyclic-fold and a hetcroclinic collision occurring simultaneously; see £|3j 

SnicHom A SNIC bifurcation and a homoclinic collision occurring simultaneously; see ^}5| 

SnicHet A double SnicHom bifurcation occurring simultaneously in Zi symmetric systems; see 

The SnicHom bifurcation is particularly important in our problem because it separates the two 
possible scenarios upon entering the pinning region: the stable limit cycle outside may disappear 
in a homoclinic collision or a SNIC bifurcation. 



8 Fluid dynamics examples of pinning due to breaking the 
50(2) symmetry 

8.1 Pinning in small aspect ratio Taylor-Couette flow 

Experiments in small aspect-ratio Taylor-Couette flows have reported the presence of a band in 
parameter space where rotating waves become steady non-axisymmetric solutions (a pinning effect) 
via infinite-period bifurcations (Pfister et al. 1991). Previous numerical simulations, assuming 



5*0(2) symmetry of the apparatus, were unable to reproduce these observations ( Marques fc Lopez] 
2006 ) . Recent additional experiments suggest that the pinning effect is not intrinsic to the dynam- 



ics of the problem, but rather is an extrinsic response induced by the presence of imperfections 
that break the 5*0(2) symmetry of the ideal problem. Additional controlled symmetry-breaking 



perturbations were introduced into the experiment by tilting one of the endwalls ( Abshagen et al. 



2008 1. Pacheco et al. (2011) conducted direct numerical simulations of the Navier-Stokes equations 



including the tilt of one endwall by a very small angle. Those simulations agree very well with the 
experiments, and the normal form theory developed in this paper provides a theoretical framework 
for understanding the observations. A brief summary of those results follows. 

Taylor-Couette flow consists of a fluid confined in an annular region with inner radius r% and 
outer radius r D , capped by endwalls a distance h apart. The endwalls and the outer cylinder are 
stationary, and the flow is driven by the rotation of the inner cylinder at constant angular speed f2 



(see figure 22 for a schematic). The system is governed by three parameters: 



the Reynolds number Re = Qri(r — ri)/v, 
the aspect ratio r = h/(r D — n), 

the radius ratio ( n = r i/ r o, 



(41) 
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Figure 22: Schematic of the Taylor-Couette apparatus. 
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(a) 




(b) 




Figure 23: (a) Experimental regimes found in the small aspect-ratio Taylor-Couette problem, with 
a pinning region, adapted from Pfister et al. (1991). (6) Numerically computed rotating wave from 
Marques & Lopez (2006), found in regions II±. 



where v is the kinematic viscosity of the fluid. The system is invariant to arbitrary rotations about 
the axis, SO (2) symmetry, and to reflections about the mid-height, a Z 2 symmetry that commutes 
with SO(2). In both the experiments and the numerical simulations, the radius ratio was kept fixed 
at r\ — 0.5. Re and T were varied, and these correspond to the parameters /i and v in the normal 
forms studied here. 

For small i?e, below the curve PF in figure [23 |( a), the flow is steady, axisymmetric and reflection 
symmetric, consisting of two Taylor vortices (Marques & Lopez 2006). The Zi reflection symmetry 
is broken in a pitchfork bifurcation along the curve PF, and a pair of steady axisymmetric one- 
vortex states that have a jet of angular momentum emerging from the inner cylinder boundary 
layer near one or other of the endwalls is born. Both are stable, and which is realized depends 
on initial conditions. The only symmetry of these symmetrically-related solutions is SO{2). The 
inset in figure [22] shows the azimuthal velocity associated with the state with the jet near the top. 
These steady axisymmetric one-vortex states are stable in region I. There are other flow states 
that are stable in this same region. For example, above the dashed curve SN in figure [23|a), 
the two-Taylor-vortex state becomes stable and coexists with the one-vortex states. However, the 
two-vortex and the one-vortex states are well separated in phase space and the experiments and 
numerics we describe below are focused on the one-vortex state. On increasing Re, the one-vortex 
state suffers a Hopf bifurcation that breaks the SO(2) symmetry and a rotating wave state emerges 
with azimuthal wave number m = 2. Figure [23^ 6) shows an isosurface of axial angular momentum, 
illustrating the three-dimensional structure of the rotating wave. For slight variations in aspect 
ratio, the rotating wave may precess either prograde (in region II + above H + ) or retrograde (in 
region II_ above H_ ) with the inner cylinder, and in between a pinning region III is observed. This 
is observed even with a nominally perfect experimental system, i.e. with the SO(2) symmetry to 
within the tolerances in constructing the apparatus. The Hopf bifurcation is supercritical around 
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Figure 24: Bifurcation diagrams for the one-cell state from (a) the experimental results of Abshagen 
et al. ( 2008[) with the natural imperfections of their system, and (b) the numerical results of Pacheco 
et al. (20111 with a tilt of 0.1° on the upper lid. The dotted curve in both is the numerically 



determined Hopf curve with zero tilt. 



the region where the precession frequency changes sign. However, for smaller aspect ratios the Hopf 
bifurcation becomes subcritical at the Bautin point Ba in figure [23{ a). 

Various different experiments in this regime have been conducted in the nominally perfect sys- 
tem, as well as with a small tilt of an endwall (Pfister et al. 1988| [1991 1992 Abshagen et al. 



20081. Figure 24 ' a) shows a bifurcation diagram from the laboratory experiments of Abshagen 



et al. (2008). These experiments show that without an imposed tilt, the natural imperfections of 
the system produce a measurable pinning region, and that the additional tilting of one endwall 
increases the extent of the pinning region. Tilt angles of the order of 0.1° are necessary for the 
tilt to dominate over the natural imperfections. Figure |24[6) shows a bifurcation diagram from 



the numerical results of Pacheco et al. (2011 ), including a tilt of one of the endwalls of about 0.1°, 
showing very good agreement with the experimental results. Included in figures 24 [a) and (b) is 
the numerically computed bifurcation diagram in the perfect system, shown as dotted curves. The 
effects of imperfections are seen to be only important in the parameter range where the Hopf fre- 
quency is close to zero and a pinning region appears. It is bounded by infinite-period bifurcations 
of limit cycles. The correspondence between these results and the normal form theory described 
in this paper is excellent, strongly suggesting that the general remarks on pinning extracted from 
the analysis of the five particular cases are indeed realized both experimentally and numerically. 



These two studies (Abshagen et al. . 2008 Pacheco et al. 2011 ) are the only cases we know of where 



quantitative data about the pinning region are available. Yet, even in these cases the dynamics 
close to the intersection of the Hopf curve with the pinning region, that according to our analysis 
should include complicated bifurcational processes, has not been explored either numerically or 
experimentally. This is a very interesting problem that deserves further exploration. 



8.2 Pinning in rotating Rayleigh-Benard convection 

Up to now, we have considered the zero-frequency Hopf problem in the context of a supercritical 
Hopf bifurcation. However, in the Taylor-Couette example discussed in the previous section, the 
zero frequency occurs quite close to a Bautin bifurcation, at which the Hopf bifurcation switches 
from being supercritical to subcritical, and a natural question is what are the consequences of 
the zero-frequency occurring on a subcritical Hopf bifurcation. The normal form theory for the 
behavior local to the Hopf bifurcation carries over by changing the direction of time and the sign 
of the parameters /x and v as discussed before, but then both the limit cycle and the pinned state 
are unstable and not observable in a physical experiment or direct numerical simulation. The limit 
cycle becomes observable as it undergoes a saddle-node of limit cycles (a cyclic fold) bifurcation at 
the fold associated with the Bautin bifurcation (see a schematic in figure [25]a,) , and we expect that 
the pinned state does likewise with a saddle-node of fixed points bifurcation along the same fold. 
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Figure 25: (a) Schematic of the Bautin bifurcation, including the path L of frequency zero bifurcated 
states in the 5*0(2) perfect system, shown in a two parameter space (/i, v) with A a global measure 
of the solution. Also shown is the projection of the saddle-node surface on parameter space, (b) 
Schematic of the rotating convection apparatus, with the streamlines of the basic state shown in 
the inset. 



In this subsection, we identify a rotating convection problem where precisely this occurs ( Marques 
et al.\ 2007 Lopez & Marques 2009), and conduct new numerical simulations by introducing an 



SO{2) symmetry-breaking bifurcation that produces a pinning region on the upper branch of the 
subcritical Hopf bifurcation. 

The rotating convection problem consists of the flow in a circular cylinder of radius ro and height 
h, rotating at a constant rate w rad/s. The cold top and hot bottom endwalls are maintained at 
constant temperatures Tq =p 0.5AT, where Tq is the mean temperature and AT is the temperature 



difference between the endwalls. The sidewall has zero heat flux. Figure 25 'b) shows a schematic 
of the flow configuration. 

Using the Boussinesq approximation that all fluid properties are constant except for the density 
in the gravitational and centrifugal buoyancy terms, and using h as the length scale, h 2 /k as the 
time scale, and AT as the temperature scale, the governing equations written in the rotating frame 
of reference are: 

. „ ^ t „ „ crFrRa , . 

(d t +u- V)w = -Vp + <tV 2 m + aRaQz + 2uQu x z (6 - z)r, (42) 

7 

{dt + u- v)e = ui + v 2 e, v-u = o, (43) 

where u = (u,v,w) is the velocity field in cylindrical coordinates (r, 6, z) in the rotating frame, p 
is the kinematic pressure (including gravitational and centrifugal contributions), z the unit vector 
in the vertical direction z, and r is the radial vector in cylindrical coordinates. Instead of the 
non-dimensional temperature T, we have used the temperature deviation with respect to the 
conductive profile, T = Tq/ AT — z + 0, as is customary in many thermal convection studies. 
There are five non-dimensional independent parameters: 

the Rayleigh number Ra = agh 3 AT / {kv) , 

the Froude number Fr = ui 2 ro/g, 

the Coriolis number Q = u>h 2 /v, (44) 

the Prandtl number a = v/k, 

the aspect ratio 7 = ro/h, 

where a is the coefficient of volume expansion, g is the gravitational acceleration, n is the thermal 
diffusivity, and v is the kinematic viscosity. 
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Figure 26: (a) Bifurcation curves for fi = 100, a = 7 and 7 = 1, where Hzp are the segments of 
the Hopf bifurcation with negative and positive frequency, the Bautin point Ba is where the Hopf 
bifurcation switches from super- to subcritical and the cyclic-fold bifurcation curve CF emerges. 
The lines L; and L u are the loci where the rotating wave has zero frequency on the lower and 
upper branches of the cyclic fold. The rectangle around L u corresponds to figure 27 1. (b) The 
frequency along the Hopf bifurcation H T with the Bautin point and the point ZF, where the sense 
of precession changes, marked as open symbols. 
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The boundary conditions for u and O arc: 



r = 7 : 
z = ±1/2 



Q r = u = v = w = 0, 
O = u = v = w = 0. 



(45) 
(46) 



For any Fr ^ 0, the system is not invariant to the so-called Boussinesq symmetry corresponding to 
invariance to a reflection K z about the half-height z — 0, whose action is K z (u, v, w, 0,p)(r, 9, z) = 
(u,v,—w, — 0,p)(r, —z). The system is only invariant under rotations about the axis of the 
cylinder, the 5*0(2) symmetry. 

The governing equations have been solved using a second-order time-splitting method combined 
with a pseudo-spectral method for the spatial discretization, utilizing a Galerkin-Fourier expansion 
in the azimuthal coordinate 9 and Chebyshev collocation in r and z. The details are presented in 



Mercader et al. (2010). We have used n r = 36, ng — 40 and n z = 64 spectral modes in r, 9 and 
z and a time-step dt — 2 x 10~ 5 thermal time units in all computations. We have checked the 
spectral convergence of the code using the infinity norm of the spectral coefficients of the computed 
solutions. The trailing coefficients of the spectral expansions are at least five orders of magnitude 
smaller than the leading coefficients. In order to compute the zero-frequency line L in the subcritical 
region of the Bautin bifurcation, where the fixed points and limit cycles involved are unstable, we 
have used arclength continuation methods for fixed points and for rotating waves adapted to our 



spectral codes (Sanchez et al. 2002 Mercader et al. 2006). 



Figure |26[ a) shows the parameter region of interest in this convection problem. In the region of 
high Froude number (region I) we have a stable steady solution, consisting of a single axisymmetric 
convective roll where the warm fluid moves upwards close to the axis (due to the rotation of the 
container), and returns along the sidewall, as illustrated in figure 25 b); the inset shows streamlines 
of the flow in this base state, that is S'0(2)-equivariant with respect to rotations about the cylinder 
axis. This base state looses stability when the Froude number Fr decreases, in a Hopf bifurcation 
along the curves H±. The bifurcation is supercritical for Ra < 14157 and subcritical for higher 
Ra; the change from supercritical to subcritical happens at the codimension-two Bautin bifurcation 
point Ba, at (Ra, Fr) ss (14 157, 0.3684). The bifurcated limit cycle, a rotating wave with azimuthal 
wave number m — 3, is unstable, but becomes stable at the cyclic fold curve CF (a saddle- node 
bifurcation of limit cycles). This curve CF originates at the Bautin point Ba. There are other flow 
states that are stable in this same region ( Lopez & Marques . 2009 1 ; these additional states are well 



separated in phase space and the numerics we describe below are focused on the base state and the 
m = 3 bifurcated rotating wave. 



Figure 26 b) shows the computed frequency of the limit cycle along the Hopf bifurcation curve. 
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Figure 27: (a) SNIC bifurcation curves in (Ra, Fr) space bounding the pinning region for fl = 100, 
a = 7 and 7 = 1. The region shown is the rectangle in figure 26 1, and the dashed line is the 
L u curve in the perfect case, (b) The period of the rotating wave C3 as it approaches the SNIC 
bifurcations at Fr = 0.32 and other parameters as in (a). 



This frequency is negative along H_ and positive along H + , and is zero at the ZF (zero frequency) 
point. At this point we have precisely the scenario discussed in the present paper: a flow (the base 
state) with 50(2) symmetry undergoing a Hopf bifurcation that has zero frequency at that point. 
Figure [26^ a) also includes the line L in parameter space where the frequency of the bifurcated states 
is zero. This curve has been computed using continuation methods since the zero-frequency state 
is unstable in the lower part (L;) of the saddle- node CF, and therefore cannot be obtained via time 
evolution. The zero frequency state becomes stable upon crossing the saddle-node curve CF and 
moving to the upper part L u of the saddle- node CF, and becomes observable both experimentally 
and by numerical simulations advancing the Navier-Stokes equations in time. 

In order to break the 50(2) symmetry and see if a pinning region appears, an imperfection has 
been introduced, in the form of an imposed linear profile of temperature at the top lid: 

Q(r,e,z) =ercos0 at z = l/2, (47) 

where e is a measure of the symmetry breaking. This term completely breaks the rotational sym- 
metry of the governing equations, and no symmetry remains. Figure |27[ q) shows that the line 
L becomes a band of pinned solutions, steady solutions with frequency zero, as predicted by the 
normal form theory presented in this paper. We can also check the nature of the bifurcation taking 



place at the boundary of the pinning region. Figure 27 [b) shows the variation of the period of the 
limit cycle approaching the pinning region. It is an infinite period bifurcation, and the square root 
fit (shown in the figure) works better than the logarithmic fit. We estimate that the bifurcation is 
a SNIC bifurcation, as the normal form theory presented predicts it should be sufficiently far from 
the zero frequency point ZF. 



Figure 28 shows contours of the temperature at a horitzontal section at mid height (z = 0) for 
the symmetric system (e = 0) in figure [28{ a), and for the system with an imperfection e = 0.05, 
corresponding to a maximum variation of temperature of 5% AT at the top lid, in figure [28[ 6). The 



parameter values for both snapshots are Ra = 21950 and Fr = 0.32, inside the pinning region in 



figure |27[a). The pinned solution 28 [b) has broken the 50(2) symmetries, is a steady solution, and 
we can see that one of the three arms of the solution is closer to the wall than the other two. The 
attachment of the solution to the side wall, due to the imperfection at the top lid, results in the 
pinning phenomenon. 



9 Summary and conclusions 

The aim of this paper has been to provide a general dynamical systems description of the pinning 
phenomenon which is observed in systems possessing two ingredients: slowly traveling or rotating 
waves and imperfections. The description boils down to the unfolding of a Hopf bifurcation in an 
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Figure 28: Temperature contours at mid height (z = 0) for ft = 100, a = 7, 7 = 1, Ra = 21950 
and Fr = 0.32. (a) is the symmetric solution without imperfection (e = 0), and (b) is a pinned 
solution with an imperfection e = 0.05. There are 20 quadratically spaced contours in the interval 
T e [-0.31,0.31], with blue (red) for the cold (warm) fluid. 

50(2) equivariant system about the point where the Hopf frequency is zero. This turns out to be a 
very complicated problem due to the degeneracies involved, but by considering all of the low-order 
ways in which SO(2) symmetry may be broken near a zero frequency Hopf bifurcation, we can 
identify a number of general features which are common to all scenarios, and hence can be expected 
to be found in practice. These are that the curve of zero frequency splits into a region in parameter 
space of finite width that scales with the strength of the imperfection, and this region is delimited 
by SNIC bifurcations. In the very small neighborhood of the zero frequency Hopf bifurcation 
point, where the SNIC curves and the Hopf curve approach each other, the dynamics is extremely 
complicated, consisting in a multitude of codimension-two local bifurcations and global bifurcations. 
The details depend on the particulars of the imperfection, but all of these complications are very 
localized and are not resolvable in any practical sense. We provide two examples in canonical 
fluid dynamics to illustrate both the pinning phenomenon and the use of the theory to describe 
it. These are a Taylor-Couette flow in which the Hopf bifurcation is supercritical and a rotating 
Rayleigh-Benard flow where the Hopf bifurcation is subcritical. 
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A Notation and description of bifurcations 



Codimension-one bifurcations 

Name Description 

SN±,o Saddle-node (also called fold) bifurcations 

H±,o Hopf bifurcations 

PF± Pitchfork bifurcations 

CF Cyclic fold: two limit cycles are born simultaneously 

L, L;.„ Limit cycle becoming a family of fixed points 

Hom± ; o, Horn Homoclinic collision of a limit cycle with a saddle 

Heti^ Hcteroclinic collision of a limit cycle with saddles 

SNIC±,o Saddle-node appearing on a limit cycle 

Glu Gluing bifurcation - two limit cycles collide with a saddle 



Codimension-two bifurcations 


Name 


Description 


Cusp± 


Cusp bifurcations 


TB±, TB 


Takens-Bogdanov bifurcations 


dPF± 


Degenerate pitchfork - zero cubic term 


Ba 


Bautin bifurcation - degenerate Hopf with zero cubic term 


PfGl 


Simultaneous gluing Gl and pitchfork PF bifurcations 


CfHom 


Simultaneous cyclic-fold CF and homoclinic collision Horn 


CfHet± 


Simultaneous cyclic-fold CF and heteroclinic collision Horn 


SnicHom± ; o 


Simultaneous SNIC and homoclinic collision 


SnicHeti^ 


Simultaneous SNIC and heteroclinic collision 



B Symmetry breaking of 50(2) with an e term: computa- 
tions 

Fixed points 



The fixed points of the normal form ( 12 ) are given by r — (f> = 0, i.e. 

2t:t--J{) - -^--^-•ow, (48) 

resulting in the cubic equation f(p) — p 3 — 2(ap + bv)p 2 + (p 2 + v 2 )p — 1 = 0, where p = r 2 . This 
equation always has a real solution with p > 0, and in some regions in parameter space may have 
three solutions. The curve separating these behaviors is a curve of saddle-node bifurcations, where a 
couple of additional fixed points are born. This saddle-node curve is given by f(p) = f'(p) = 0; from 
these equations we can obtain (/i, v) as a function of p. In order to describe the curve it is better 
to use the rotated reference frame (u, v), where the it-axis coincides with the line L, introduced in 
(13) (see also figure [2ja)). The saddle-node curve is given by 



(u,v) = ^ 1 (i + 2pI±Ja P 1-i), p 2 e (2- 2 / 3 ,+ TO ), (49) 



2p 2 



where pi is the double root of the cubic equation f(p) = 0. The third root p is given by pop% = 1< 
If there is any point where the three roots coincide (i.e. a cusp bifurcation point, where two saddle- 
node curves meet), it must satisfy f(p) = f'(p) = f"{p) = 0. There are two such points Cusp±, 
given by p 2 — 1 and (it, u)cusp ± = (3/2, ±v3/2), dividing the saddle-node curve into three branches: 
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SN , joining Cusp + and Cusp_, and unbounded branches SN + and SN_ starting at Cusp + and 
Cusp_ respectively, and becoming asymptotic to the line L. Along SNo, P2 < 1 < Po, while along 
SN + and SN_ , po < 1 < p2 ■ At the cusp points, the three roots coincide and their common value 
is +1. 

A better parametrization of the saddle-node curve is obtained by introducing s — ± \J (4/9?, — l)/3, 
so that now s G (— oo, +oo), Cusp-t corresponds to s = ±1 and 



(u,v) 



(3(1 



Po 



[2(1 + 3s 2 )] 2 / 3 ' 
Hopf bifurcations of the fixed points 



(l + 3s 2 ) 



2/3 



P2 



l + 3s 2 \V3 



(50) 



Using Cartesian coordinates z = x + iy in ( 12 ) we obtain 



p —b 
v p 



(x 2 + y 2 ) 



ax — by\ 
bx + ay J ' 



where we have set e = 1. The Jacobian of the right-hand side of (51) is given by 

j _ (p — 3ax 2 — ay 2 + 2bxy —v + bx 2 + 3by 2 — 2axy\ 
\ v — "ibx 2 — by 2 — laxy p — ax 2 — Say 2 — 2bxy J 

The invariants of the Jacobian are given by 

T = 2(p — 2ar 2 ), D = p 2 + v 2 - 4(ap + bv)r 2 + 3r 4 . 

A Hopf bifurcation takes place iff T = and D > 0. When T = 0, 4a 2 D = 4(c 



bp? 



(51) 



(52) 



(53) 



The fixed points satisfying T — are given by f(p) — and p — 2ap, resulting in the curve T in 
parameter space 

(54) 



Aapv{av — bp) = 8a 3 — p 3 , 



that can be parametrized as 



( / x^)=a 1 /3(l- s 2)V 3 ^2,^ + 



vT 



s e (-1.+1) 



(55) 



For s — > ±1, p = and v — > ±oo, and the curve is asymptotic to the p = axis, the Hopf curve 
for e = 0. Along the T curve, the determinant D is given by 



D = a 2 / 3 (l - s 2 )- 1 / 3 ^ 2 - 1 - 2^sVl^ 2 ): 
resulting in two Hopf bifurcation curves (when D > 0): 



H_ 



s e 



(_l,-V(l-&)/2), H + : se (v/(l + &)/2,+l). 



(56) 



(57) 



The end points of these curves have T = D = 0, and are Takens-Bogdanov bifurcation points TB±. 
They are precisely on the saddle-node curve (50), where both curves are tangent. The coordinates 
of the four points Cusp± and TB± are: 



(M> ^)cusp 4 
(m,^)tb+ 



b , a \ 

a p,H — , 

V3 y/3' 



{p,v)cusp_ = - [a+-j=,b- 



TB„ 



(2a, 2b- 1) 



,1/3 ■ 



(58) 
(59) 



(2a, 26+ 1) 

(2(1 + 6)) 1/3 ' " ' (2(1-&)) J 

As TB± are on the saddle-node curve, its s± parameter, according to (|50|), can be computed. The 



result is s s = 8y/{l - <56)/3(l + 6b), with 8 = ±'. Therefore TB + G SN , closer to Cusp+ than to 
Cusp_; TB_ G SN_ if a < 60°, TB_ G SN when a > 60° and TB_ = Cusp_ for a = 60°. 
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C Symmetry breaking of 50(2) to Z^'. computations 
Fixed points 



The fixed points of the normal form ( 19 ) are given by r = 6 = 0. One solution is r = (named 



Pq). The other fixed points are the solutions of 

cos 26 — ar 2 — a. 1 , 2 \2 , , i\i -, /™\ 

s^26 = v- br 2 \ =* (M-ar 2 ) 2 + (,-6r 2 ) 2 = l, (60) 

resulting in the bi-quadratic equation r 4 — 2(a/x + bv)r 2 + fi 2 + v 2 — 1 = 0, whose solutions are 

r| = an + bv ± A, e 2i0± = (av - b[i =F iA)e iQ ° (61) 
A 2 = (a/i + bv) 2 + 1 - ^ 2 - v 2 = 1 - (av - fc^) 2 , (62) 

and every admits two solutions, differing by n (they are related by the symmetry Z 2 , z — > —z 
discussed above). Introducing a new phase ot\, 

av~bfi-iA = c iai , (63) 

where a± G [— it, 0] because A > 0, we immediately obtain: 

e 2i0± = e i( Q0 ±«i) ^ ± = ( ao ±a 1 )/2 1 (64) 

with the other solution being (oto ± ai)/2 + 7r; a.\ is a function of (/i, i/) while ao is a fixed constant. 
We have obtained two pairs of Z 2 symmetric points, P + = r + e 1 ^+ and P£ = —r + e 1 ^ + , and P_ = 
r-e^- and PI = -r_e^-. 

Hopf bifurcations of fixed points 



The Jacobian of the right-hand side of ( 26 1 is given by 



j _ At + 1 — 3ax 2 — ay 2 + 2bxy —v + bx 2 + 3by 2 — 2axy \ , g, 
I ^ — 36x 2 — fey 2 — 2axj/ fi—1 — ax 2 — 3ay 2 — 2bxy J 

The invariants of the Jacobian are the trace T, the determinant D and the discriminant Q = T 2 —4D. 
They are given by 

T = 2(/i — 2ar 2 ), (66) 

D = n 2 + v 2 - 1 - 4(a/i + bv)r 2 + 3r 4 + 2(a(a; 2 - y 2 ) - 2bxy) , (67) 

Q = A(l - v 2 + Abvr 2 + (1 - 46 2 )r 4 - 2(a(x 2 - y 2 ) - 2bxy)}. (68) 



The eigenvalues of the Jacobian matrix ( 65 ) in terms of the invariants are A± = h (T ± y/Q) . For 
example, a Hopf bifurcation takes place iff T = and Q < 0. For the fixed points P s and P*, where 
s = ±, we obtain 

T(P S ) = T(P*) = 2((b 2 - a 2 )^ - 2abv - 2asA), (69) 
D(P s ) = D(P:)=4sAr 2 , (70) 
Q(P S ) = Q(P S *) = 4(( M - 2ar 2 ) 2 - 4 S Ar 2 ). (71) 

As a result, Q(P ) = Q(P*) > and P_ and P* never experience a Hopf bifurcation. After some 
computations, T(P + ) — results in the ellipse 

/i 2 - Aabfiv + 4a 2 v 2 = 4a 2 , (72) 
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centered at the origin, contained between the straight lines SN + and SN_ and passing through the 
points (fj,,v) — (0, ±1), the ends of the horizontal diameter of the circle /j, 2 + v 2 = 1. This ellipse 
is tangent to SN + and SN_ at the points (/i, v) = ±(26, (b 2 — a 2 ) /a). The condition Q < is only 
satisfied on the elliptic arc from (fj,, v) — (0,-1) to (26, (b 2 — a 2 )/a) with fx > 0; along this arc 
P + and undergo a Hopf bifurcation. We have assumed that a and b are both positive. The 
properties of the ellipse are: 

major semiaxis (1 — ^-)/U = 2abv, length 2a/y / IL, (73) 

minor semiaxis (1 — = — 2abfi, length 2a/\/£+, (74) 



where 2£± = 1 + 4a 2 ± y/l + 8a 2 . The eccentricity e is given by 



2 , l+4a 2 

= l + -f= - ■ 75 

e 2 Vl + 8a 2 



Codimension-two bifurcations of fixed points 

The Jacobian evaluated at the three points TB + , TB_ and TB is: 



J(TB + )=(| (76) 
J(TB_)=(\ (77) 
J(TB) = ( * ( 1 + ;V a V (78) 



(b-l)/a -1 

The three matrices have double-zero eigenvalues and are of rank one, so the three of them correspond 
to Takens-Bogdanov bifurcations. The state that bifurcates at the TB point is P+, without any 
symmetry, so that it is an ordinary Takens-Bogdanov bifurcation, although the Z2 symmetric state 
P* also bifurcates at the same point in parameter space (but removed in phase space) at another 
ordinary Takens-Bogdanov bifurcation. The state that bifurcates at the TB± points is P$. This 
state is Z2 symmetric, and so these are Takens-Bogdanov bifurcations with Z2 symmetry. 
The Jacobian evaluated at the two points dPF± is 

,<dPF + ) = (Y -Z b ). W = C_V <™> 

The corresponding eigenvalues arc A + = —26 and A_ = for dPF + , and A + = 26 and A_ = for 
dPF_. Both points are pitchfork bifurcations, and in order to determine if they arc degenerate, their 
normal form needs to be computed in order to verify that the cubic term is zero. However, since in 
a degenerate pitchfork bifurcation a curve of saddle-node bifurcations emerges that is tangent to the 
pitchfork bifurcation curve, from figure [7] it is immediate apparent that both dPF-t are degenerate 
pitchfork bifurcations. 

D Symmetry breaking of 50(2) with quadratic terms: com- 
putations 



The three cases (|32j), (35) and (30) can be dealt with by considering the normal form 

f = r(u — ar 2 ) + r 2 cos md>, 

■ , - (80) 

© = v — br + r sin mm. 



where m = 1 for the cz case (^5.3), m — — 1 for the ezz case ([5.2) and m — —3 for the ez case 



(q5T|). The fixed points, other than the trivial solution Pq (r = 0), in the three cases are given by 
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the biquadratic equation r 4 — 2(a/i + bv + l/2)r 2 + /i 2 + v 2 = 0, with solutions P± 



<\ = afi + bv + 1/2 ± (an + bv + 1/4 - (ai/ - bfi) 



2\l/2 



(81) 



The phases <j> of the P± fixed points can be recovered from 

cos m(j> — ar — /i/r, sinm0 = fer — vjr. (82) 

For rn = ±1 the solution is unique; for m — 3 the solutions come in triples, differing by 27r/m. 
It is convenient to use the phase space coordinates adapted to line L, introduced in ( 13 ) (see also 



figure |2[a)). In terms of these coordinates, r\ 



1/2 ± y/u + 1/4 — v 2 , and the fixed points P± 
exist only in the interior of the parabola u — v 2 — 1/4, whose axis is the line L. On the parabola, 
these points are born in saddle-node bifurcations. In order to explore additional bifurcations of 
these points, we compute the Jacobian matrix of the normal form (80), 

J 



ji — 'iar 2 + 2r cos mc/ 
—2br + siamd) 



-mr sinmtj 
mr cos md> 



whose trace and determinant, for the P± points, are easily computed: 
T(P±) = (m - l)ar 2 ± - (to + 1)//, 

D(P ± ) = -2m^u + 1/4 - v 2 (^J(u + 1/2) 2 - u 2 - v 2 ± (u + 1/2)) . 



(83) 



(84) 



Therefore signD(P-i-) = =p sign to, and for to > (to < 0) only P + (P_) may undergo a Hopf 
bifurcation. 



Th 



e ez case 



({5.1). Here to = — 3 and T = 2/i — Aar 2 . P_ is a saddle, but P + undergoes a 



Hopf bifurcation when T = 0. The condition T = for P + gives 

v/a/i + + 1/4 - (av - bfi) 2 = — -((a 2 - b 2 )fi + 2abfi + a) 

2a V / 



> 0. 



(85) 



By squaring and simplifying, we obtain the ellipse (bfj, — 2av) 2 + (a/i — l) 2 = 1 which is tangent to 
the line fj, — at the origin, with its center at (/x, v) = (2a, b)/(2a 2 ), and whose elements are: 



major semiaxis parallel to (1 — = 2abv, 
minor semiaxis parallel to (1 — lJ)v = —2abfi, 



length 

length l/^/l^, 



(86) 
(87) 



where 2£± — 1 + 4a 2 ± vl + 8a 2 . This ellipse has much in common with the one found in the ez 
case, and the eccentricity e is given by the same expression (75 1 . For > 7r/6, the ellipse is located 
in the interior of the parabola of saddle-nodes, for ao — tt/6 it becomes tangent to the parabola at 
a single point, and for ao < tt/6 it becomes tangent at the two points 



H = - (l - 2a 2 - sb\J 1 - 4a 2 ) , 
a v ' 

v= ^v 7 ! - 4a?(by/l - 4a 2 - s(l - 2a 2 )), 



s = ±1. 



These are the points TB S in figure 15 b). Only the points on the elliptic arc H joining these two 
points satisfy (85), and along this arc P + undergoes a Hopf bifurcation. 

The ezz case ({ 5.2). Here rn = —1 and T(P±) = —2ar 2 < 0, so there are no Hopf bifurcations. 
Moreover, D(P + ) > 0, so it is always stable and L*(P_) < 0, so it is a saddle. The only exception 
is when r = 0, and this only happens at /i = v = 0, the degenerate high-codimension point at the 
origin. 

The ez 2 case ({5.3). Here to = 1, and T = — 2/x is zero on the line ji = inside the parabola. 
On this line Ho, P+ undergoes a Hopf bifurcation, and the points of contact with the parabola have 
D = T = so they are Takens-Bogdanov bifurcations (see figure [T7]b) . The Hopf and Takens- 



Bogdanov bifurcations are degenerate, as will be discussed in { D.l 
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D.l A degenerate Takens-Bogdanov bifurcation 



In the ez case, numerical simulations of the normal form ( 34 ) show that the Hopf bifurcation Hq and 



the Takens-Bogdanov points TB-t are degenerate. This can also be found by direct computation. 
Let us work out the details for the TB_ point. 

The coordinates of TB_ in parameter space are (fi, v) = (0, —0.5/(1 + 6)) = (0, 1/4 cos 2 (a/2)), 
where P± are born in a saddle-node bifurcation, and the fixed points are given by 



-ia/2 



2(1 + 6) 4cos 2 (a/2) 



z± 



2(1 + 6) 4cos 2 (a/2) 



(89) 



In order to obtain the normal form corresponding to the Takens-Bogdanov point, a translation of 
the origin plus a convenient rescaling of z and time is made: 



t = 4rcos 2 (a/2), C = 2(z - z±) cos(a/2), p, + w = 4(/i + iv) cos 2 (a/2). 



Substituting in (34) results in 



( = iC + ie" 2iQ C + e ia/2 C 2 + 2e" 3iQ/2 |C| 2 -ie^ iQ CIC| 2 - 
In order to obtain the normal form, we introduce the real variables [x\,y\) 

C=( yi +2i Xl )e- ia , 

so that the linear part of the ODE is transformed into Jordan form, and we obtain 



(90) 



(91) 



(92) 



cos(a/2) 



2x12/1 
4x 2 + 3y 2 



sin(a/2) 



-2x\ ~ 3j/?/2 



(Ax 2 + y 2 ) 



x sin a + (y/2) cos a 



(93) 



y sin a — 2x cos a 

Now we can reduce the quadratic and cubic terms to normal form by an appropriate near-identity 



quadratic transformation (x\,yi) — ¥ (x2,y2)- Knobloch (1986) gives explicitely the normal form 
coefficients up to and including third order, in terms of the coefficients of the original ODE (in the 
form 93); a nice summary is also given in Wiggins (2003 §19.9). Using this explicit transformation, 
we obtain 

±2 = 2/2 

y 2 = 4x| cos(a/2) + 16cc| cos 2 (a/2) + 0(4) 



(94) 



and the xiyi term in the normal form of the Takens-Bogdanov bifurcation is missing, resulting in 
a degenerate case, the so-called cusp case, of codimension three. The unfolding of this degenerate 



case has been analyzed in detail in Dumortier et al. { 1987). Note that the ODE (94) is Hamiltonian 



at least up to order three, which helps to explain the continuous family of periodic orbits obtained 



in the interior of the homoclinic loop in figure 17 6). 
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